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Abstract 

Searches for supersymmetry are among the most exciting physics goals at Run II of the Tevatron. 
In particular, in supersymmetric models with light charginos, neutralinos and sleptons, associated 
chargino-neutralino production can potentially be observed as multi-lepton events with missing en- 
ergy. We discuss how, in the generic TeV-scale MSSM, the prospects for these chargino-neutralino 
searches are impacted by cosmological considerations, namely the neutralino relic abundance and di- 
rect detection limits. We also discuss what an observation of chargino-neutralino production at the 
Tevatron would imply for the prospects of future direct dark matter searches without assuming specific 
patterns of supersymmetry breaking. 
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I. INTRODUCTION 



A consensus has formed within the astrophysics community in support of the conclusion that 
the majority of our universe's mass takes the form of cold, collisionless dark matter Despite 
the very large body of evidence in favor of dark matter's existence, the nature of this elusive 
substance remains unknown. Of the many dark matter candidates to have been proposed, one 
of the most compelling and most often studied is the lightest neutralino in i?-parity conserving 
models of supersymmetry [2|. 

Among the most prominent missions of the Tevatron's Run II are its searches for supersym- 
metry. Results from Tevatron searches for squarks and gluinos [sj, neutralinos and charginos [3], 
stops and sbottoms [5i], and the Higgs bosons of the Minimal Supersymmetric Standard Model 
(MSSM) [ol have each recently been published. While no evidence for supersymmetry has yet 
been found, in many cases these results represent the strongest limits to date. Although the 
Tevatron is not well suited to directly place limits on the properties of the lightest neutralino, 
the results of these other searches can have considerable implications for the nature of this dark 
matter agent. 

One of the prime channels for observing supersymmetry at the Tevatron is associated 
neutralino-chargino production 0, [8|. These particles can decay to the lightest neutralino 
and leptons through the exchange of either sleptons or gauge bosons, resulting in events fea- 
turing three leptons and missing energy. The results of these searches are somewhat model 
dependent, but the current results from CDF and DO can be used to exclude charginos as 
heavy as approximately 150 GeV in some models, well beyond LEP's chargino mass limit of 
104 GeV. By the end of Run II, the Tevatron is expected to exclude selected models with 
charginos not far below 200 GeV. 

Neutralino dark matter can be detected through its elastic scattering with nuclei. Exper- 
imental efforts designed to observe such events are known as direct detection. The prospects 
for this class of techniques depends on the composition of the lightest neutralino, as well as 
on the masses and couplings of the exchanged squarks and Higgs bosons. Generally speaking, 
information from collider searches for supersymmetry, whether detections or constraints, can 
be used to better estimate the prospects for the detection of neutralino dark matter. The 
relationship between Tevatron and LHC searches for heavy MSSM Higgs bosons and direct 
searches for neutralino dark matter has been studied in detail elsewhere |9|. Here, we return to 
this theme, but focussing on searches for trilepton events from associated neutralino-chargino 
production at the Tevatron (see also Ref. [lo|). 

In the past, it has been possible to link the Tevatron's trilepton signature to other signa- 
tures for new physics, for example the decay Bg fifi [11]. Such links rely for example on 
a correlation between light sleptons and small values of tan (3 for the chargino and neutralino 
decays on one hand and the pseudoscalar Higgs boson mass and large values of tan (3 in flavor 
physics on the other. Naively, similar correlations should be present when the dark matter 
candidates annihilate mainly through an s-channel Higgs boson and the trilepton signature 
requires relatively light supersymmetric scalars. Moreover, one could imagine correlations be- 
tween these signals in the co-annihilation region, if the lightest slepton is mass degenerate with 
the lightest neutralino, limiting the visibility of the trilepton channel. However, these fairly 
obvious correlations rely on a series of assumptions. First, the different MSSM scalar masses 
have to be correlated. Secondly, the light scalar masses should in some way be linked to the 
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FIG. 1: The cross section for associated chargino-neutralino production at Tevatron Run II, as a 
function of the lightest chargino mass, for various choices of tan/3 (3 and 60, in the left and right 



frames, respectively) and [Ij]. In these Figures we use 2Mi 
At = Ab = Ar =0. 
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lightest neutralino mass and to the mass difference between the hght chargino and neutrahno. 
Last but not least, the dark matter particle should annihilate dominantly through one channel. 
The aim of this analysis is to determine how much of a correlation between dark matter and 
Tevatron searches survives if we assume only a TeV-scale MSSM spectrum with no specific 
patterns of supersymmetry breaking. 

This article is structured as follows: In Section [TTl we discuss the searches for associated 
chargino-neutralino production at the Tevatron. In Section llllt we turn our attention to the 
thermal relic abundance of neutralinos, focussing on those models within the reach of the Teva- 
tron and the correlation between Tevatron measurements and the neutralino's relic density. 
In Section HVt we discuss direct detection prospects for such models and the correlations be- 
tween those and Tevatron observations. Finally, we summarize our results and conclusions in 
Section |Vl 



II. NEUTRALINO CHARGINO SEARCHES AT THE TEVATRON 

In many supersymmetric models, associated chargino-neutralino production can occur with 
a cross section on the order of a picobarn at Run II of the Tevatron (1.96 TeV center-of-mass 
collisions). These particles can each subsequently decay to the lightest neutralino and leptons 
(Xi^ — >■ Xi ^^^1 X2 ~^ Xi l^l^)^ either through the exchange of charged sleptons or gauge bosons. 
This can lead to distinctive trilepton plus missing energy events which, in some supersymmetric 
models, could be identified over Standard Model backgrounds. 

In order for SUSY-trilepton events to be extracted at the Tevatron, however, the underlying 
supersymmetric model must possess a number of rather specific features. In particular, the Xi^ 
and X2 must both be light. In Fig. [T] we plot the associated chargino-neutralino production 
cross section as a function of the lightest chargino mass for various values of tan (3 and /i. The 
cross section drops rapidly for heavy chargino/neutralino masses. Additionally, in order to be 
identified at the Tevatron, Xi^ and X2 decays must each occur with large branching fractions to 
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FIG. 2: The branching fractions of xi ^-nd Xi decays to final states with charged leptons, as a function 
of the Ughtest stau mass and for tan /3 = 3, 60. The three fines denote fi =200 GeV (dashed), 500 GeV 
(dotted) and 1 TeV (sofid). We have also used 2Mi = M2 =140 GeV, rUg = nif^ = tua =500 GeV and 
At = Ah = At =0. In order for the combination of X2xf to decay mostly to final states with three 
charged leptons, the lightest stau (possibly along with other charged sleptons) must be rather light. 

charged leptons, which means that the supersymmetric mass spectrum is arranged such that 
chargino and neutrahno decay primarily to charged sleptons rather than to (off-shell) gauge 
bosons or squarks, each of which lead to significant branching fractions to jets. Furthermore, 
WZ production leads to a Standard Model background of trileptons plus missing energy from 
which any SUSY-trileptons must be separated. To accomplish this, the analyses of CDF and 
DO each include kinematic cuts on observables like mu, designed to remove backgrounds. They 
reduce the efficiency for supersymmetric events with charginos and/or neutralinos decaying 
through gauge bosons essentially to zero. 

To ensure large branching fractions for charginos and neutralinos through slepton exchange, 
the lighter sleptons must be quite light. To avoid large neutralino or chargino branching frac- 
tions to neutrinos, the sneutrino masses (along with the left handed charged sleptons) must be 
somewhat heavier. Unless we want to break the SU{2) symmetry between charged slepton and 
sneutrino masses, this means the lightest charged slepton should be dominantly right handed, 
independent of unification assumptions. In Fig. [2l we plot the branching fractions of charginos 
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FIG. 3: The current limits on associated neutralino-chargino production from CDF and DO searches 
for SUSY-trilepton events. See text for details. 

and neutralinos to trileptons (including electrons, muons and taus), as a function of the lightest 
stau mass, for various choices of /i and tan p. 

Limits from the CDF and DO collaborations have been placed on the combined cross section 
for associated neutralino-chargino production and branching fractions to three leptons. In 
particular, DO has published results for their search for events with three leptons (at least two 
of which are electrons or muons) plus missing energy using the first 320 pb~^ of data from 
Run II [l^. They find a rate consistent with the predictions of the Standard Model, and use 
this to place constraints on supersymmetry. CDF has published the results of their search for 
events with two like-sign leptons (electrons or muons) and missing energy using 1 fb~^ of data 
from Run II In this analysis, 13 events were observed, a slight excess compared to the 
7.8 predicted by the Standard Model (corresponding to a chance probability of 7%). More 
recently, CDF has published the results of their combined search for associated neutralino- 



chargino production. These findings are consistent with Standard Model expectations [15|. In 
addition to these published results, a number of preliminary results from CDF [16| and DO 17 
searches for trilepton plus missing energy events have been reported. 

In Fig. [3|, we show the current limits from CDF and DO in this channel. The limits from 
CDF are shown for two scenarios, labelled "mSUGRA" and "no-mixing". In the mSUGRA 
scenario, the masses of the staus are determined within the context of the mSUGRA model, 
which leads to the lightest stau being considerably less massive than the other sleptons and, in 
turn, to large branching fractions for chargino and neutralino decays to taus. In the no-mixing 
scenario, decays to taus, muons and electrons are approximately equally common. As taus 
are more difficult to identify than other leptons, the CDF limit in the mSUGRA scenario is 
considerably weaker than in the no-mixing case. Also shown is the DO limit for the no-mixing 
scenario. It is more stringent than the limit from CDF, in part, because DO's result is slightly 
stronger than expected. By the end of Run II, the limits from each of these experiments are 
expected to improve by a factor of approximately 5 to 10. 
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III. THERMAL ABUNDANCE OF NEUTRALINOS WITHIN TEVATRON REACH 



In i?-parity conserving models in which the hghtest neutrahno is the hghtest supersymmetric 
particle (LSP), such particles fall out of thermal equilibrium when the rate of Hubble expansion 
begins to dominate over their annihilation rate. The resulting density of neutralino dark matter 
in the universe today is related to its annihilation cross section: 

^yOh^ ^ — -, (1) 

where (av) is the thermally averaged neutralino-neutralino annihilation cross section, g* is 
the number of relativistic degrees of freedom available at the temperature of freeze-out, and 
xp = m^o/Tp, where Tp is the temperature of freeze-out. For neutralinos (and other species of 
electroweak scale WIMPs), xp falls in the range of 20-30. The thermally averaged annihilation 
cross section can be written as {av) ~ a -f- Sb/xp, where a and b are terms in the expansion 
av = a + bv^ + 'd{v*). 

The neutralino annihilation cross section depends on the details of the supersymmetric 
model, including the composition of the LSP and the masses and mixings of the exchanged 
sparticles and Higgs bosons. The four neutralinos of the MSSM are mixtures of the superpart- 
ners of the photon, Z and neutral Higgs bosons. The neutralino mass matrix is diagonalized 
into mass eigenstates by a unitary rotation N*M^oN~^. Hence, we can describe the lightest 
neutralino as a mixture of gauginos and higgsinos: 

X? = A^n^ + A^isiy^ ^ j^^^fj^ ^ (2) 

Although no accelerator bounds have been placed on the mass of the lightest neutralino 
directly [iSj, LEP II has placed a lower limit of 104 GeV on the mass of the lightest chargino, 
which is in turn related (at tree level) to M2, tan/? and /i. 



1 



M2P + I^P + 2m^ - J(|M2|2 + |/i|2 + 2ml.y - A\fiM2 - m\ sin2/5|^ 



1/2 



(3) 



The LEP II bound, therefore, leads to a constraint of IM2I, > 104 GeV. Since we are 
interested in the case in which the xf ci-nd X2 within the reach of the Tevatron, and yet 
significantly heavier than the lightest neutralino, we are forced to consider values of Mi smaller 
than M2 and If Mi is considerably smaller than M2 and \fi\, the lightest neutralino will be 
largely bino-like, with a small higgsino admixture: 

|A^i3p ^ m| sin^6'i4' sin^/3 
liViiP \l^ 

The mass of the lightest neutralino in this scenario is approximately given by 

m| sin^ 6w{Mi + jj sin 2(3) 
IJ,^ - Ml 



/200GeVy 






\Nu\^ < \Ni3\ 







kr^-— — ^ --o.oi( ,,, ) , liv^r < liViar . (4) 



"^X? ^ ^1 -2 ■ (5) 



In most supersymmetric models within the reach of trilepton searches at the Tevatron, 
the lightest neutralino typically annihilates somewhat inefficiently and thus is expected to be 
produced in the early universe with a thermal abundance in excess of the measured dark matter 
density. There are a number of possible exceptions to this conclusion, however. In particular: 



6 




FIG. 4: The thermal neutrahno rehc abundance as a function of its mass for various values of the 
slepton masses. The slepton masses neglecting mixing are (from top to bottom) 1 TeV (black), 
300 GeV (red), 200 GeV (blue), 140 GeV (green) and 120 GeV (magenta). Other parameters are 
niA = 'mq=l TeV and tan/3=10. M2 is either 2Mi or the lowest value consistent with the LEP 
chargino bound, whichever is greater. The trilinear couplings At and Ai, are selected to maximize 
the light Higgs mass. In the left (right) frame, as set to 5Mi (lOMi). The horizontal dashed lines 
denote the dark matter abundance measured by WMAP [19]. The two dips correspond to the Z and 
light Higgs resonances. 



- If the lightest neutralino is within a few GeV of the Z 01 h resonances (2m^o ^ mz^h), 
then annihilations through these channels can be very efficient, especially if the neutralino 
has a sizable higgsino fraction {i.e. moderate to small values of |/i|). For example, the 
cross section for Z-mediated neutralino annihilation scales simply as the square of the 
difference of the two higgsino fractions, (lA'^iap — |A'^i4p)^. Its effect can be seen in Fig. |H 

- Light sleptons, which are required in models within the reach of trilepton searches at the 
Tevatron, can also lead to efficient neutralino annihilation. In the extreme case, the light- 
est stau can be quasi-degenerate with the lightest neutralino, leading to highly efficient 
coannihilations. The effect of sleptons in the neutralino relic abundance calculation can 
be seen in Fig. |H In this figure, we show the relic density as a function of the LSP mass, 
for various values (1000, 300, 200, 140 and 120 GeV) of the slepton masses. ^ 

- If the currently (largely) unconstrained pseudoscalar Higgs boson is light enough 
and its couplings are large (large tan/3 and/or small |/i|) then it will efficiently mediate 
neutralino annihilations. When not near the A^-resonance, the cross section to down-type 
fermions through pseudoscalar Higgs exchange is proportional to Mf ta.ia^ (3 m'j/\fi\m\. 
This contribution is most significant in the case of a mixed gaugino-higgsino with a light 
pseudoscalar Higgs and large tan /3. 

From Fig. HJ it is obvious that light neutralinos will be overproduced in the early universe unless 



By "slepton mass" or mf, we refer to a common mass for the selectrons, smuons and staus before off-diagonal 
terms in the mass matrices are accounted for. This quantity approximately corresponds to the selectron and 
smuon masses. The staus, in contrast, will depart somewhat from this value, m? m~ =F to,-(^t — /itan/3). 
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FIG. 5: Left panels: the thermal neutralino relic abundance as a function of the neutralino mass, in 
models within the reach of Tevatron trilepton searches. Dark points denote models which have already 
been excluded by the Tevatron trilepton searches, whereas the lighter points are within reach with 8 
fb^^ of integrated luminosity. In the left frame, all of the models within the Tevatron reach are shown. 
In the center frame, we omit models with efficient Z, h or ^''-mediated dark matter annihilation (see 
text for details). In the right frame, we show the cross section times branching ratio for trilepton 
production at the tevatron as a function of the thermal relic abundance of neutralinos. In this frame 
we only show models with LSP masses within the range of 70zbl GeV. 



the sleptons are light, the lightest neutralino's mass is within a few GeV of the Z or h resonances 
or pseudoscalar Higgs exchange provides a significant contribution to the annihilation cross 
section. 

We demonstrate this further in Fig. El where we compare the relic abundances found in 
various models within the Tevatron reach. In this parameter scan, we vary the masses Mi, M2, 
rrii, rriq, \fi\ and rriA up to 1 TeV. Values of tan/5 within the range of 1 to 60 are considered. For 
simplicity, we assume the gluino mass to be M3 ^ 3.7M2. All models shown in Fig. [5] satisfy 
all collider constraints on the chargino, slepton, squark and Higgs masses. The relic abundance 
we compute using DarkSUSY [2Q] . 

In the left frame of Fig. [5] we show all models found to be within the Tevatron reach. In 
the second frame we omit models with neutralinos annihilating through either a. Z,h reso- 
nance or via an s-channel diagram. To quantify this selection we remove parameter points 
in which the lightest neutralino mass is within 7 GeV of the Z or h resonances, or with 
(tan /3/10)^/[(m^/l TeV)^ (1/^1/1 TeV)] > 1. These cuts overwhelmingly remove models with 
low relic densities, thus demonstrating that neutralino annihilation through Z or h resonances 
or through A^-exchange are essentially required if dark matter is to avoid being overproduced 
in models within the reach of the Tevatron. Note however, that independent of the LSP mass 
well always find models which produce the correct relic density. 

In the right frame of Fig. [5] we fix the LSP mass to 70 ± 1 GeV and show the correlation 
between the trilepton cross section times branching ratio versus the relic density. The mass of 
the produced neutralino and chargino is free. The fact that the majority of points tend towards 
overdosing the universe corresponds to a bias in the entire data sample, also seen in the left 
two panels of the same figure. We checked that independent of the LSP mass chosen there is 
indeed no visible correlation between the relic density and the Tevatron trilepton cross section 
in the MSSM. 
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IV. DIRECT DETECTION OF NEUTRALINOS IN TEVATRON REACH 



Experiments such as XENON HH, CDMS and many others [23| have over the last several 
years placed increasingly stringent limits on the elastic scattering cross section of WIMPs with 
nuclei. The neutralino's elastic scattering cross section with nuclei is given by 



4m ^0*^' 



TT [771^0 + TTIt) 



[Zfp + {A-Z)U 



(6) 



where rriT is the target nuclei's mass, and Z and A are the atomic number and atomic mass of 
the nucleus, fp and /„ are the neutralino's couplings to protons and neutrons, given by: 



g=u,d,s 



q=c,b,t 



(7) 



where Ug are the neutralino-quark couplings and f^J ^ 0.020 ± 0.004, f^^' ^ 0.026 ± 0.005 



/^^^ ^ 0.118 ± 0.062, f^^' ^ 0.014 ± 0.003, ^ 0.036 ± 0.008 and f^^' ^ 0.118 ± 0.062 [24 

The first term in the above equation corresponds to interactions with the quarks in the target 
nuclei, whereas the second term denotes interactions with the gluons in the target through a 



,(n) 



,(n) 



(n) 



dp) 



quark/squark loop diagram. is given by 1 

/S ^ 0.83. 

The neutralino-quark coupling is given by [2 



ip) 



0.84, and analogously. 
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(8) 



(9) 



In these expressions, i = 1,2 denote up and down- type quarks, respectively. mii,m2i denote 
the squark mass eigenvalues and 1] is the matrix which diagonalizes the squark mass matrices. 
i/i, Tsi and denote hypercharge, isospin and electric charge of the quarks. For scattering off 
of up-type quarks 6i = Ni3,62 = Ni4,B = sin/5, C = sinQ;,i5 = cos a, whereas for down-type 
quarks 6i = Ni4,62 = —Ni3,B = cos/3,C = cosa,D = — sin a. a is the mixing angle in the 
Higgs sector. 
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FIG. 6: The neutralino's elastic scattering cross section with nucleons, as a function of its mass, in 
model which are within the reach of trilepton searches at the Tevatron with 8 fb^^ luminosity. In 
both frames, the dark points represent models which are predicted to generate a thermal density of 
dark matter within the range measured by WMAP [19|]. In the left (right) frame, the lighter points 
represent models which predict a larger (smaller) dark matter density than is measured. 

In Fig. [6l we plot the neutralino-nucleon elastic scattering cross section found in those 
models which are within the 8 fb~^ reach of the Tevatron's trilepton search. In each frame, 
the dark points correspond to models which predict a thermal abundance of neutralino dark 
matter within the range measured by WMAP ^]. The lighter points represent models with too 
much (left) or too little (right) dark matter relative to the measured abundance. As expected 
from the different composition patterns of the neutralinos, we find a very large range of cross 
sections, varying from about 10~^ to 10~^^ pb. Models which predict an abundance of dark 
matter below the measured value tend to have somewhat larger elastic scattering cross sections 
with nuclei. This is due to the coupling to heavy MSSM Higgs bosons which can both mediate 
neutralino annihilation and elastic scattering processes. Most models with the observed relic 
density fall in the upper portion of the elastic cross section range. In particular, the majority 
of them are within the reach of CDMS's current run (labelled CDMS 2007). This is not the 
case for a typical scan over the entire MSSM parameter space, which consists mostly of models 
beyond the Tevatron's reach (see, for example, Ref. (26|). 

The reasons for the favorable direct detection prospects among models within the reach of 
the Tevatron are somewhat subtle. In the case of a neutralino annihilating in the early universe 
primarily through a. Z or h resonance, little can be said regarding the prospects for direct 
detection. Furthermore, in models which annihilate largely through slepton exchange in the 
early universe (or through coannihilations with sleptons), the elastic scattering cross section is 
likely to be suppressed. In many of the models within the reach of the Tevatron, however, the 
neutralino annihilation cross section is dominated by pseudoscalar Higgs exchange. In these 
models, which feature moderate to large values of tan P and somewhat light pseudoscalar Higgs 
masses, the elastic scattering cross section is typically dominated by the exchange of the heavy 
scalar Higgs, H, with strange and bottom quarks, leading to a neutralino-nucleon cross section 
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of: 

The similarity rriH ^ ttla is of course not an artefact of a SUSY-breaking assumption, but a 
generic feature of the two-Higgs-doublet model. Because in a large fraction of these models 
the combination of tan^ /3/m^|/i| is large in order to generate an acceptable relic abundance, 
the direct detection rates also have a tendency to be larger compared to those found in a more 
general sample of supersymmetric models. 

V. OUTLOOK 

In this article, we have studied the cosmological implications of supersymmetric models 
within the reach of searches for associated neutralino-chargino production at Run II of the 
Tevatron. We have analyzed how results from this Tevatron search channel might impact the 
prospects for direct searches for neutralino dark matter. Although there is not a particularly di- 
rect or obvious connection between these two experimental programs, it is important to consider 
how to exploit the interplay between collider and astrophysical searches for super symmetry. 

Supersymmetric models whose rate of trilepton events from associated chargino-neutralino 
production is within the reach of the Tevatron have some rather peculiar features. In particular, 
they contain light neutralinos which either annihilate through a Z or h resonance, through 
pseudoscalar Higgs exchange, or via very light sleptons. Therefore, for models in which the 
lightest neutralino is not within a few GeV of mzf^ or m/i/2, the heavy Higgs bosons A,H 
tend to be light and values of tan/3 are typically moderate to large. These neutralinos will 
also show a non-negligible higgsino fraction. These features lead to larger than average elastic 
scattering cross sections with nuclei (which is dominated by H exchange), and high rates in 
underground direct dark matter experiments. This means that if the Tevatron detects trilepton 
events from associated neutralino-chargino production, the near future prospects for the direct 
detection of neutralino dark matter are expected to be promising. 

The absense of more distinct parameter correlations means that the collider and the cosmo- 
logical analyses of the MSSM neutralino and chargino sector probe different properties of the 
supersymmetric spectrum. This is different from the case of, for example, gravity-mediated 
SUSY breaking. Looking at the TeV-scale MSSM this implies that the information gained in 
dark matter searches is largely orthogonal to the information which could be obtained from 
collider searches. Only by combining many sets of information from many different experimen- 
tal channels will it become possible to construct with confidence a consistent picture of the 
TeV-scale Lagrangian [23]. 
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